ON THE NON-HOMOGENEOUS NAVIER-STOKES SYSTEM WITH 
NAVIER FRICTION BOUNDARY CONDITIONS 



LUCAS C. F. FERREIRA, GABRIELA PLANAS, AND ELDER J. VILLAMIZAR-ROA 

Abstract. We address the issue of existence of weak solutions for the non-homogeneous 
Navier-Stokes system with Navier friction boundary conditions allowing the presence of 
vacuum zones and assuming rough conditions on the data. We also study the conver- 
gence, as the viscosity goes to zero, of weak solutions for the non-homogeneous Navier- 
Stokes system with Navier friction boundary conditions to the strong solution of the 
Euler equations with variable density, provided that the initial data converge in to a 
smooth enough limit. 



1. Introduction 

We are concerned with the incompressible Navier-Stokes model with variable density 
in a bounded domain. The governing equations are given by the following system 

dt{pu) + div (puu) — uAu + Vvr = pf in Q, 

div u = in Q, (1.1) 
dtp + div (pu) = in Q. 

Here, Q = x (0, T), where is a bounded domain of M'^ with smooth boundary dO., and 
T > 0. The unknowns are the velocity field u, the density p, and the pressure vr of the 
fluid. The parameter v > is the viscosity coefficient of the fiuid and / is a given vector 
field driving the motion. 

We supplement the system (|l.ip with initial and Navier friction boundary conditions 

n • n = on S, 
[D{u)n + au]^^^ = on S, , 

p{o) = po in n, ^'■^> 

{pu){0) = vq in ri, 

where S = dQ x (0,T), n is the exterior normal vector to di^, po > denotes the initial 
density and vq has to be at least such that V(){x) = whenever po{x) = 0. Moreover, 
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D{u) = ^{diUj + djUi)i<ij<n denotes the deformation tensor, [■]tan is the tangential com- 
ponent of a vector on dVl, and puu = p{u® u). The constant a > stands for the friction 
coefficient which measures the tendency of the fluid to slip on the boundary. 

The goal of this paper is to study the convergence of solutions for ()l.ip - ()1.2p . as the 
viscosity goes to zero, toward the solution of the Euler equations with variable density. 
Formally, when we drop the viscous term (i.e., taking = 0) system (jl.ip - ()1.2p degenerates 
into the non-homogeneous Euler equations 



dt{pu) + div (puu) + Vvr = 


Pf 


in 


Q 


div u 


= 


in 


Q 


dtp + div (pu) 


= 


in 


Q 


u ■ n 


= 


on 


s 


piO) = 


= Po 


in 


n 


(H(o) = 




in 





(1.3) 



We aim at giving a justification of this formal procedure. 

The issue of the vanishing viscosity limit or inviscid limit for the incompressible homo- 
geneous Navier-Stokes equations is a classical problem in fluid mechanics. In the whole 
space and periodic cases, the inviscid limit was performed by several authors, see e.g. 
[H El [T8| I35j . In the case where there exist physical boundaries, the problem of conver- 
gence leads to the formation of a boundary layer if one supplements the Navier-Stokes 
equations with no-slip boundary conditions (which are the most often prescribed ones). 
This happens because there is a discrepancy between the no-slip boundary conditions for 
the Navier-Stokes equations and the tangential boundary conditions for the Euler equa- 
tions. 

There is no consensus on the boundary conditions to be prescribed for the Navier-Stokes 
equations, except for impermeable boundary which corresponds to the condition (11.2P i . 
Navier |28] claimed that the tangential component of the viscous stress at the boundary 
should be proportional to the tangential velocity, leading to the boundary condition (I1.2p 9. 
Conditions (jl.2P i- (|1.2p 9 are called Navier friction boundary conditions, or simply Navier 
boundary conditions (another names have been used as well). These conditions were also 
derived by Maxwell |26j from the kinetic theory of gases and rigorously justified as a 
homogenization of the no-slip condition on a rough boundary (see |17j). 

Recently, the inviscid limit for the Navier-Stokes equations with Navier boundary con- 
ditions was established, for which the reader is referred to [6l 113 1 l23 l 125]. The situation in 
this case is thus very different from the case of no-slip boundary conditions and requires 
distinct involved arguments. 

Concerning the non-homogeneous incompressible Navier-Stokes equations, it is worth- 
while to remark that there exists a considerable number of papers devoted to their mathe- 
matical analysis, especially in the case where the equations are complemented with Dirich- 
let boundary conditions. Those results can be classified in two classes: on the one hand, 
there are existence results when the initial density is assumed to be positive and so there is 
no vacuum initially; and on the other hand, the case where the initial-vacuum is allowed. 
The first case has been addressed by several authors, see e.g. [H [71 [191 [211 [291 [311 [3l [15] , and 
references therein. In order to avoid vacuum, the basic assumption in the above-quoted 
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works is 

< Co < inf po{x) < po{x) < suppo(2;), 

and so, in particular po, p have a positive lower bound. In the second one, when the 
initial- vacuum is allowed (po ^ 0), the problem ()l.ip - ()1.2p is more difficult to handling. 
Indeed, comparing with the first case, fewer results are available in the literature related 
to the existence of weak solutions (see [32l [33l [201 122] )• In particular, to the best of our 
knowledge, existence of weak solutions has not been still treated for the non-homogeneous 
incompressible Navier-Stokes system with Navier boundary conditions. Let us mention 
the work [16] where strong solutions to the system with slip boundary conditions are 
considered, however vacuum zones are not admitted. Thus, our first goal will be show the 
existence of weak solutions for ()l.ip - (jl.2p . allowing vacuum and assuming rough conditions 
on the field / and initial data po and vq. 

The vanishing viscosity limit for the non-homogeneous incompressible Navier-Stokes 
equations in the whole space, or with periodic boundary conditions, was addressed in 
|14 1 I15 1 [8]. as long as no vacuum states occur. They proved the convergence of local strong 
solutions in Hilbert and Sobolev spaces. One of the difficulties in this case is to show that 
the time existence is independent of the viscosity. As well as in the homogeneous case, it is 
expected that the issue of the inviscid limit in bounded domains presents boundary-layer 
phenomenon when one considers no-slip boundary condition. 

Our second goal in this paper is to show that a weak solution of the non-homogenous 
incompressible Navier-Stokes equations (II. ip with Navier boundary conditions (II. 2p con- 
verges in the energy space toward the strong solution of the non-homogeneous incompress- 
ible Euler equations ()1.3p in the inviscid limit. This extends some earlier results obtained 
by [U [23l [m [30] in the homogeneous case. The strategy is to compare the smooth solution 
of the Euler equations and a weak solution of the Navier-Stokes equations, which is the 
leading idea of the proof of weak-strong uniqueness for the homogeneous incompressible 
Navier-Stokes equations. In our case, this approach arises suitably in view of the fact that 
only strong solutions are known to exist for the non-homogeneous Euler equations (11. Sp 
(see e.g. |34]). In fact, to best of our knowledge, there is no theory of weak solutions for 

This paper is structured as follows. In the next section, we establish a result of existence 
of weak solutions with finite energy for the problem (|l.ip - (jl.2p (see Theorem 12. 3p . The 
precise definition of weak solution with finite energy is given in Definition 12.11 below. In 
Section [3l we address the vanishing viscosity limit for the non- homogeneous Navier-Stokes 
equations with Navier friction boundary conditions (see Theorem 13. 2p . To this end, we 
first recall a result of local strong solution for the Euler equations with variable density 
(see Theorem 13. ip . Finally, we proceed with the proof of the inviscid limit. 

We finish this section by establishing some notations used throughout this manuscript. 
We denote by P(il) and P'(il) the space of functions of class C°°{Q) with compact support, 
and the space of distributions on il, respectively. We use standard notations for Lebesgue 
and Sobolev spaces. We denote by || • ||p the norm in U'(i}). Otherwise, the norm will be 
specified. For a Banach space X, we indicate by (•, •)x',x the duality product between X' 
(the dual space of X) and X. As usual, we will use the same notation for vector valued 
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and scalar valued spaces. There will be no danger of confusion since the difference will 
be clear in the context. Also, we denote by H^{^) the subspace of H^{Q,) of divergence 
free vector fields tangent to the boundary. Finally, C^([0,T];X) represents the space of 
functions it : [0, T] — )• X which are continuous with respect to the weak topology. 



2. Weak solutions for the non-homogeneous Navier-Stokes system 



In this section we study the existence of global weak solutions for the non-homogeneous 
incompressible Navier-Stokes with Navier boundary conditions. We assume that the initial 



2p_ 



I 1 2 

density po belongs to L^(f]), 6 < p < oo, allowing to vanish, vq € Lp+^(Q) with G 

L^{Q,) and the external force / G L^{0,T; L~ (Q.)). The data vq and correspond 

formally to initial value for pu and p\u\ , respectively. In the case p = oo, L''-exponents 
depending on p should be understood in the natural way, that is, as the limit when p — )• oo. 

For instance, Lp+i and Lp-i become L in that case. 

Now we introduce the definition of weak solution with finite energy for the system 

(inD-dLi. 

Definition 2.1. A weak solution for ()l.ip - ()1.2p is a pair of functions {u, p) verifying the 
following items: 

i) u e L^{0,T-H^{n)), p G C{[0,T];W~^'P{n)) n L°°{0,T;LP{n)), p>0 a.e. in 

Q, pu G L°°(0,T;L^(O)), ^u G L^{0,T; L^{n)), puu G L\0,T; L^{n)), such 
that equation dtp+ div (pu) = is satisfied in T)'{Q) and the momentum equation 
^TT^i is verified in the following sense: 



/ I pu - dtif + 2az/ / 
Jo Jn Jo 



u ■ (f + 2u 



an 



T 



D{u) : D{^) 

■T 







puu • Vif, 
pf-ip+f vo^iO), (2.1) 



for all (f e C^{[0,T];Hl{n)), ip{T,x)=0 a.e. in Q. 
a) The initial data (|1.2p q is verified in the following sense: 



Hi) The following energy inequality 



Po^l^dx, V V G w^'^'in). 



^\\,/^u{t)\\l+2ua 



\u\ 



2v 



Jan 



\Du\ 



< - 



vo 



Jn 



pf-u (2.2) 



holds for a.e. t G (0,T). 



Let us make some commentaries about the previous definition. We first notice that the 
divergence- free boundary condition of the velocity field and the boundary condition (jl.2p i 
are given by the choice of the space H^{Q,). The weak formulation (|2.ip also contains 
the boundary condition (jl.2p 9 in the sense that if u is more regular, say H^{i^) H H'^{VL), 
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it can be recovered. Indeed, first let us just recall, for the readers convenience, that the 
formulation (j2.ip comes from following integration by parts: 



Lemma 2.2. (^|13j ) Let f and g be smooth vector fields such that g is divergence free and 
tangent to the boundary. Then 

- [ Af-g = 2 [ D{f) : D{g) - 2 / [D{f)n]tan ■ 9- 
Jn Jn Jan 

Now, assuming that u is more regular and using the previous lemma, from ()2.ip . we 
obtain 

[ [ [D{u)n + au]-(p = 0, 
Jo Jan 

for any test function satisfying ip ■ n = on dil.; consequently, [D{u)n + au]tan = on 

an. _ 

We also note that by taking test functions in 'P([0, T) x in the form iph = '4^{x)0h{z), 
where V S T){Q) with divergente free, and 9^ G 'D([0,T)) such that 9h{z) = 1 for z < t 
and Oh{z) = for z > t + h, and taking the limit as /i — )• 0, we obtain an equivalent weak 
formulation for the momentum equation 

{pu){t)ij + 2au [ [ u-i) + 2u [ [ D{u):D{ip)- [ [ puu-Vip 



Jan Jo Jn Jo Jn 

Jn Jn 

2p 

from which we deduce that pu S C^([0, T]; Lp+i (17)). Hence, the initial data (ll.2P /i is 
verified in the following sense: {pu){t) converges weakly to vq as t — t- O"*". 

The result of existence of weak solutions with finite energy is the following. 

Theorem 2.3. Let6<p<oo, f e L'^{0,T; {Q)), po G LP{n), po > 0, a.e. inU, and 

m 



Vq G Lp+i(r2), -1^^ G L^{VL). There exists a weak solution {p,u) of problem (|l.ip - ()1.2p 



po 



the sense of Definition \2.1\ 

Remark 2.1. Let us observe that the density p belongs to C([0, T]; L^(il)). In fact, noting 
that u G L2(0,r;F^(f])) C L^{Q,T]L^{VL)) and I < p' < 2 < Q, the desired claim follows 
from standard approximation and regularization arguments for transport equations (see 
e.g. [9]j. On the other hand, when po{x) > cq > 0, the weak solution of (|l.ip - (|1.2p given 
by Theorem {2^ also verifies u G L°°{0,T; L'^{il.)). 

To prove the existence of a weak solution with finite energy, we first introduce a reg- 
ularized problem, depending on a small positive parameter e, which is constructed by 
a regularization of the continuity and momentum equations, as well as, a regularization 
of the data. More explicitly, fixed e > 0, we consider the following regularized problem 
related to (|l.ip - (|1.2p : Find (ne,/Oe) solution of system 
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dt{peUe) + div [p^UeUe) - v/S.Ue + VvTe = Pef + ^{Ap^)Ue ill Q, 

div = in Q, 
dtPe + div (peUe) = cAp^ ill Q, 
tig • n = on S, 
[D{u,)n + an4^„ = on S, (2.3) 
dpe 
dn 

/5e(0) = Po,e in f^, 
(/9en,)(0) = i;o,e in 0, 

where /Oo,e, similarly to [HI p. 149] for 6 < p < cxo, is such that /3o,e G C^'^(il), r G (0, 1) 
with 



on S, 



5n 



on < e < po,e(2;), x G 
|{x G r2 : pQ^e{x) < Po{x)}\ as e 0, 
Po,e — ^ Po in -^^''(^^) as e — )■ 0, 6 < p < oo, 
po,e ^ weakly- * in L°°(0) as e — )• 0, p = oo, 

and the initial linear momentum vq ^ is defined as 



(2.4) 



vo if PoA^) ^ Po(a:), 
if poA^) < Poix)- 



(2.5) 



We introduce the concept of weak-strong solution to the previous regularized system. 



Definition 2.4. Let 6 < p < oo. A weak-strong solution of (12. 3p is a pair of functions 
{u„p,) satisfying G L^{0,T; HA^)), Pe G C([0, T]; P^-1'P(J7)) nL-(0, T; LP(0)), > 
a.e. in Q, p,u^ G L~(0, T; (17)), peii^Ue G L1(0, T; 1.2(17)), Vp^ G ^^(0, T; L2(17)), 

i) equation p.SP '^ /lo/ds a.e. m Q. T/ie boundary condition (12. Sp ^ /loWs a.e. on S, 

and i/ie initial condition (j2.3P 7 /lo/ffs a.e. in $7, 
iij ^/le momentum equation (I2.3p i is verified in the following sense 



tp + 2v 



an 

T 



D{u,) : D{^) 
T 



- [ [ p,u,u,-V^ + ^ [ [ Vp,V{u,-^)= [ [ pj-ip+ [ voM^): (2-6) 
Jo Jn ^ Jo Jn Jo Jn Jn 

for if G C^([0, T] X $7), with div ip = 0, ip ■ n = on Ti, and ip{T,x) = in 

In order to show the existence of weak-strong solution to the regularized problem we 
will consider the Galerkin approximations for the momentum equation and then will use 
a limiting procedure. The existence of solutions for this approximate problem will be 
obtained by linearization and the Schauder fixed point theorem. We can now state the 
result of existence weak-strong solution for the regularized system (|2.3p . 



Proposition 2.5. Let po,e be as in ()2.4p and p,uo, f as in Theorem \2.3[ Then, there exists 
{ue, Pe) a weak-strong solution of (]2.3p . in the sense of Definition \2.4\ Moreover, (u^, Pe) 

6p 

verifies p,u, £ L^{0,T; (ft)), p, G L<^{0,T;W^'<^{n)) and dtPe G L^(0, T; for 
some ^> > |- 

Proof. We split the proof into five steps. 

Step 1: Approximate problem. Let be a smooth basis of H^{Vt), or- 

thonormal in L^(il) and let 3^™ = span{w^, Let us consider a sequence {/'"jmeN 

in C([0,r];L^^(O)) such that ^ / in ^^(0, T; L^^(fi)). 

For each m E N consider the problem of finding G C([0, T]; (7^(0)) and u™ G 
Ci([0,r];3^'") satisfying 



9f/9™ + div (p™n™) = eAp™ in Q, 

— — = on S, 
an 

^"^(0) = po,e in 1^, 



(2.7) 



{dt{p'^u^) ■ V - p'^u^u'^ ■ Vv - p'^f"' ■ V + li^Du"" : Dv} dx 



+2av 



an 



[ vp"v(n'" ■v),yvey"', 

Jn 



where ulf^ G 3^™ is uniquely determined by 



Po,eUo (pax 



VQe(t>dx for ah G 3^"". 



(2.9) 



For the sake of simplicity, we shall omit the dependence of {u"^,p"^) on the parameter 
e. Observe that it™ is well-defined because the matrix with coefficients po^^w^ w"^ dx is 
invertible. Moreover, by using (|2.9p . the definition of vq^^ (see (|2.5p ) and the assumption 

|2 

G L^{^), the following estimate holds true: 



f vo,eWo*(ix = f po^e\uQ^\'^dx < f ^ '^'^^ dx 
n ' Jn ' Jn Po,e 



< 



\vo\ 
n Po 



-dx < C. 



(2.10) 



Step 2: Existence of solutions to the approximate problem. Fixed m G N, the 
existence of approximate solutions (u™,/)™) of (j2.7p - (j2.8p is proved by linearization and 
the Schauder fixed point theorem. In fact, fixed m G N and given vu G C([0, T]; 3^™), the 
following lemma gives the existence of G C{[0,T]; C^{n)) such that 

dtp"" + div (p^'w) = eAp"" in Q, 

^=OonS, (2.11) 
on 

p'-(O) = po,e in 0. 
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Lemma 2.6. ( \12\ Lemma 3.1]) Let w G C([0,T];3^™) be a given vector field. Suppose 
that po,e £ C'^'^{^), r G (0,1), irdx^o, Po,e{x) > and satisfies the compatibility condition 
9po,E _ Q Then problem ()2.1ip possesses a unique classical solution 



dn 



= p"^{w) GW = {p^G C{[0,T];C'^^in)), dtp"" G C7([0,T];C7"'^(J7))}. 



Moreover, the mapping w — t- p^{w) maps bounded sets in C([0, T]; 3^"*) into bounded sets 
in W and it is continuous with values in C^([0, T] x Q). 
Finally, as div{w) = 0, it holds 



inf po,e{x) < p"'{x,t) < suppo^eix), t G [0,r], X e^. 



(2.12) 



Posteriorly, given w G C([0,r];3^™) and G C{[0,T];C^{n)), we solve the following 
linear problem: 

Find u'^{x,t) = YlTLiMtWix) satisfying 

{dt{p"V) ■ V - p"'u'''w ■ Vv - p^'T ■ V + 2vDu'^ : Dv] dx 

+2au [ u'" -v + ^ [ Vp™V('u™ ■v) = 0,yve 3^™, 
Jan 2 Jq 

which is equivalent to 

{p'^idtu'^ + {w ■ V)^™ - D ■ V + 2i/Z)'u™ : Dv} dx 

+2av [ u'"-v + ^ [ V/)™v(u'" ■v) = o,yve y"", 

J do. 2 

n'^(O) = u"^. 

Notice that (|2.13p leads us to a linear system of ordinary differential equations for 



(2.13) 



E + E 4 WV-i + «f (*) = 0' "1 (0, T), 1 < i < m, 

{V'i(0)}?Li = components of Ug^, 



(2.14) 



where 



,2 



p'^w^w'dx G c^([o,r]). 



{(/)'"'(/; • Vw^) ■ w' + 2vDw^ : Dw'}dx + 2au / ■ 

an 



+ i f Vp<^V{uj^ . w')dx G C\[0,T]), 
2 Jn 

p'^rw'dx G cn[o,r]). 
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Since > e, it holds 

i,j "^^ i=l i=l 

It follows that the matrix A = {a\j}i^j is symmetric and positive definite, and in particular 
A is invertible. From the classical theory of ordinary differential equations, system (j2.14p 
has a unique solution {V'jjjLi e (Cn[0,r]))"; then, the solvability of the system (I2J3D 
is guaranteed. 

Multiplying equation ()2.1ip i by ^lu^p, integrating on Q and adding the result to (j2.13p 
with V = u"^{t), after some calculations, we obtain 



ld_ 

Thus, the Holder inequality implies 
1 d 



'^u'^{t)\\l + 2u\\Du"'it)\\l + 2ua [ {u'^f = [ p"^ ^ ■ vl^ 

Jan Jn 



Using a generalized Gronwall lemma ( \32\ Lemma 5]), from last inequality we can obtain 

Jo 

From the previous inequality, (j2.12p and (j2.10p . we conclude that is bounded in 
C([0,T];L2(0)) (independent of m and w). Thus, as u"'{x,t) = TljLi'^At)'^^^) and 
E^lil^iWP = \W'^\\li then {V'j} is bounded in C([0,r]) which implies that is 
bounded (independent of w) in C([0, T]; 3^™). Moreover, if w is bounded in C([0, T]; J^'"), 
from (j2.14p and the symmetry oi A = {ajj}, the set {dtipj ■ 1 < j < m} is bounded in 
C([0,T]) which implies that dtu"" is bounded in C([0,r];y'"). Thus, we conclude that u'" 
is bounded in CH[0,r];3^™) provided w is bounded in C([0, T]; 3^"^). 

Given w G C([0, T]; 3^"^), from Lemma [Ml there exists a unique p'" € C{[0,T];C'^'''{n)) 
solution of (j2.1ip . Knowing w and there exists a unique solution n™ of ()2.13p in 
C^{[0,T]]y"^) which is bounded (independent of m) in C^([0, T]; 3^™), provided w is 
bounded in C([0, T]; 3^"). Thus, there are Mi, Ma > such that \\u"'\\cmo,T];y^) < if 
ll^llc([o,r];y") < -^1- Denote by Bi the closed ball in C([0, T]; 3^"^) of radio Mi, and 
the closed ball in C^([0, T]; 3^"^) of radio M2. Then, the mapping 

T : 5i ^ ^2 

is continuous. The Arzela-Ascoli theorem implies that B2 C C([0,T];3^™) compactly and 
therefore the mapping 7~ is continuous and compact from Bi into Bi. Then, the Schauder 
fixed point theorem implies the existence of a fixed point for a given T > 0. Taking p™ 
the corresponding solution of (j2.1ip . we obtain the existence of an approximate solution 
(u™,p"^) of (fZTD -ra. 
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Step 3: Estimates for (n™, p™). We will obtain several estimates for the approximate 
solution (n"^,/?™) which are independent of m and, in general, are also independent of 
e > 0. In the sequel, C will denote a constant independent of m and e that may change 
from an estimate to another. If C depends on e we shall indicate this dependence by C(e). 

The first estimate comes from Lemma 12.61 



0< inf po,e{x) < p"'{x,t). 

Next, we multiply (I2.7p i by where p > 2, integrate by parts and use the 

incompressibility of the flow u"^ to obtain 

1 ^ M ■m / .\ Mn. / -.\ f \ \ t — I nrt i i nrt i P ^ i 



\p-^{t)\\l + e{p-l) [ \\Vpn\pn'^\' = 0, 



p dt 
then 

\\p"'\\L'^(0,T;LP{n)) < ll/50,e||LP(n) < C, (2-15) 

because (j2.4p ci and (j2.4P /i . Also, we have that 



+ 26 f \\Vp"^{s)\\lds = \\po4l (2.16) 
Jo 

so, it follows that 

iV^Vp"", V^Ap"^) is uniformly bounded in L'^{0,T; L^{n) x W-^'^{n)). (2.17) 

Integrating equation (|2.7P i with respect to the space variable and using that ^^|e = 0, 
yield the total mass conservation 

/ /3™(t)dx= / po^.dx, te[0,T]. (2.18) 
Jq Jq 



Notice that equation (|2.8p is equivalently to 

/ {p"'{dtu'^ + (n"^ • V)n"^ - /"^) • v + 2iyDu"' : Dv} dx 
Jn 

+ 2au f u'^-v-- [ V/9"^V(n"^ . t;) = G 3^™. (2.19) 
Jan 2 Jq 

Multiplying equation ()2.7P i by i|u™p, integrating on Q and adding the result with (j2.19p 
for V = u^{t), we find after some calculations 

-■^\\y/^u"'{t)g + 2u\\Du"'{t)\\l + 2ua [ jn^p = /" p™/"^ • n'". (2.20) 
2 dt Jqq Jq 

Thus, the Holder inequality implies 

~\\^u"^ml < \\Vp^u"^m2\\VF'rm2- 

Applying a generalized Gronwall lemma ( \32\ Lemma 5]), it follows that 



Jo 
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thenceforth, by applying again Holder's inequality, and using ()2.10p and (j2.15p . we obtain 



< 



II I 11^*" 11-*-/^ ll^"l|l 

PO,eUo\\2 + \\P IIl^(o,T;Lp(Q))II-/ II 



Li(0,T;LP^(n)) 



2p 



because (/'") is bounded in L^{0,T; L~ (Q)). 

On the other hand, by integrating (j2.20p on (0, t) and proceeding similarly as before, 
we arrive at 



+ 2iJ f \\Du"^\\l + 2ija f [ |u™|2 
Jo Jo Jan 



< 2llVmI^o"lli + 



1 1/2 



2p 



< 2\\VP07eU0\\2 + C\\vF^u'^\\L-°{0,T;LHn))\\p'^\\L^(OT^^^^ 

Therefore, (L>ti™) is bounded in L^{0,T; L^{n)). To estimate (n^) in L^{0,T; H^{n)) we 
shall apply the following generalized Korn inequality (see |12^ Th. 10.17]) 

\\vfHHn)<C{K,M,p){\\D{v)\\l + \\Rv\\l) 

< C{K,M,p) {\\D{v)\\l + \\R\\i\\Rv%) , (2.21) 

for V G H^{Q,) and any function R > such that < M < Rdx, \\R\\p < for some 
p > 1. Without loss of generality we can assume that po ^ and po ^ 0. It follows from 
(|2T8]1 and />o,e ^ Po in LP{^), 6 < p < oo, that 



M 



- IIpoIIi < / p"''dx = / /Oo.eC^a; < 2 || 

Polli 7 for small e > 0, 

2 Jn 



so we can take R = p™ in (|2.2ip and use previous estimates in order to infer that (um) 
is bounded in L'^{0,T;H^{il.)) C L'^{0,T;L^{Q)), independently of m and e. In the case 
p = oo, we also can choose the same M since po,e, Po > and the weak-* convergence 
/?o,e ^ Po in L'^ implies that po,tdx — p^dx. 
Next, by applying Holder's inequality, we estimate 



/ m "III 

'P^'li 12p 

L2{0,T;L^T2F(n)) 



< 



m||l/2 ll„,"i|| 

llL°°(0,T;LP(n))ll^ lll,2(0,T;L6(n)) 



< C. 



Since p > 6, we have that > 4, hence we can use interpolation to obtain 



6+2p 



/7^„,m||(p-6)/(4p-6)|| /-^^ m||3p/(4p-6) 

F II2 IIVf " II i2p 

6+2p 



Thus, taking ? = -^^5-^, it follows that ? > 2 and 

|l /777r„,m|K ^ |L/;^,,m||2(P-6)/3P|| /-^ m||2 
llY/3 u 1I4 ^ \\v P II2 II VP ^ II i2p ■ 

6+2p 

Consequently, (y^p^u™) is bounded in L'^(0, T; L^(r2)), which implies that 



^m„ m„,m| 



IP U U ||L^/2(o,T;L2(n)) 



< c. 



(2.22) 
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By using that (/o'"), [^sJJ^u!^) and (u"*) are bounded in L°°(0, T; L'^{0,T; L^{n)) 

and L^(0, T; L^(il)) respectively, we deduce the following bounds for {p"^u"^), 

"^'"^'""L^(0,T;L^(n)) - ll/''"llt(0,T;L.(O))ll^^-'"llL-(0,T;L^(^^)) < (2-23) 

II „m„,m II ^ II „m II II„,'ti|| ^ /'00/l^ 

"L2(0,T;L^^(n)) llL-(0,T;Lf{Q))l|w ||L2{0,r;L6(n)) < (^-^4) 

By interpolation it follows that 

\\P U \\2 < IIP U II 2p \\P U II 6p , 

L°°(0,T;LpTT(Q)) ^ 



where < 6^ = |^ < 1. Thus, 



2{2p - 3) 



\\p'^u"'\\Li{o,T;L^n)) < C", with g = "'"'^^ > 6. 

By using the maximal regularity for parabolic equations, we obtain that (p™) is bounded 
in Li{0,T;H^{n)), which implies that (Vp™) is bounded in L'i{0,T; L^{n)). This bound 
together the estimate for u"^ in L'^(0, T; L^(f])) allow to apply the classical L^— theory of 
parabolic equations (see [IH Th. 10.22]) to conclude that 

\\dtp"'\\Le{0,T;Le{n)) + l|P™llLe(0,T;W2,e(n)) < C(e), with Q = > |. (2.25) 

Finally, we are going to estimate the derivative in time of p"^ and p^vJ^ . To this end, 
consider the equation dtp^ = — div {p^u^) + eA/?™ in Q, and notice that W~^''^{Q) C 

VF"^'i^(0), then, by using (I2.17P and (I2.23p . we conclude that 

Wdtp'^W 1^ <C(||p™n'"|| 2, +V^||V^Ap"'||i2(oT-H/-i.2(c))) <C, 

(2.26) 

for small e > 0. 

From the momentum equation in (12. Sp together with (j2.22p . the Holder inequality and 
Sobolev imbedding, for all v G D(i7) and s > 3, we have 

4 / pVt; < ||p'"n"n™||2||Vt;||2 + 2z.||Z)^z'"||2||D^;||2 + ||p"||p||/"1l^lbll^ 

+ 2ai.||n-||^.(an)||HlL2(sn) + ^llVp'^lblln'^llellVHIa + ^llVp-lbUVn-lblblloo 

p-i 

Therefore, by using (j2.17p . we obtain 

d 



dt 



[ p'^u"'v <hm\\v\\w^,srn),yveV{n), 

Jn 

for some hm G L^(0,T), which implies that 

mp'^unh^o,T;W-^,s'^n))<C. (2.27) 
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Step 4: Convergence properties. From the uniform estimates obtained in the 
previous step, we will deduce some convergences for the approximate solution. We first 

observe that from (l2l6D - (f2T7l) . IHMh . and since H^{n) L'^{n) ^ by 

applying Lemma 4 of [32], we conclude that 



(p™) is relatively compact in L^(0, T; L^(r2)). 

Similarly, as (p™) is bounded in L°°(0, T; andLP{n) W-^'P{n) ^ W~'^"^{n), 

we have that 

(p™) is relatively compact in C{[0,T];W-^'P{Q)). 

Moreover, as s > 3 > ^1^, it holds Lh^(O) W~^'^{n) ^ W^^'"' {ft); so, from 

I^TM^i and (1X771) . we get 

[p^u"^) is relatively compact in L^{0,T;W 'p+(^{Q)). 

Thus, in view of the uniform bounds obtained in the previous step, we have that the 
sequence (n™, /9™n"^n'") converges (up to subsequences) to some (ue, Pej Cij ^2), 
as m — )• 00, in the following sense: 

p"^ Pe strongly in L'^{0, T; and a.e. in Q, (2.28) 

p"" Pe weakly- * in L~(0, T; LP{n)), (2.29) 

p™ ^ p, strongly in C([0, T]; ^^"^'^(0)), (2.30) 

-u"^ ^ -u, weakly in ^^(0, T; //^(J^)), (2.31) 

weakly in L'^{0,T; {n)), 
pm^m^^^ J weakly- * in L°°(0,r;LFK(O)), (2.32) 
{ strongly in L2(0,r;VF >+6(J7)), 



pm^m^m _^ ^^^j^jy L^/^ (0, T; (f]) ) . (2.33) 

Next, we identify the limits and .^2- To this end, notice that the product mapping 

from H^{n) X W-^^P{n) to W 'p+6(0) is continuous (see [32l Lemma 3]). Therefore, 
(p:30]l and (|2:3TD lead us to 

/9™n"* ^ weakly in L\0,T;W~^"^ (n)). (2.34) 

Convergence (|2.34|) together with (|2.32|) implies that = PeU^. In analogous way, from 
(f2:3T]l and ^3^, we obtain that 



pm^m^m ^ ^^^^^^ ^^^j^^y -^^ (0, T; P^"^' ^fe (0)) . 



The uniqueness of the limit in the sense of distributions implies that .^2 = PtUeUe- 
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Step 5: Passing to the limit as m — )• oo. Finally, we will prove the existence of a 
weak-strong solution to problem (j2.3p by passing to the limit as m — )• oo in the approximate 
problem ([TTD - OrS]) . 

Since — )■ p^j p^u"^ — )• p^Ue and A/>™ — )• Ap^ in T>'{Q), we can pass to the limit in 
(j2.7P i in the sense of distributions. Therefore 

dtPe + div (peUe) = cAp^ in V^Q). 

In view of the estimate (I2.25p . and since pe inherits regularity from p™, we have in par- 
ticular that 

dtPe + div (peUe) = eApe a.e. in Q. 

From ([2301) we have that p™(0) pe{0) in VF-^'P(ri). But, the initial condition (H^ls 
says that />"^(0) = /9o,e G C^'^(fl), so that Pe{0) = Po,e- It is not difficult to see that the 
boundary condition is also satisfied. 

In order to pass to the limit in (j2.8P i as m — )• oo, we notice that by using the energy 
identity (j2.16p for p™ and the corresponding one for p^, we can prove that Vp™ converges 
strongly to Vp*^ in L^(0, T; L^(ri)), (see [10] for a similar argument). Hence, from the 
convergences obtained in step 4, we can classically pass to the limit in (j2.8P i as m — t- oo, 
and obtain that ()2.6p holds true. The proof is complete. □ 



Proof Theorem [231 

Proof. Let {pe,Ue) be the weak-strong solution of ()2.3p given by Proposition 12.51 As we 
mentioned in step 3 of the proof of Proposition 12.51 most of the obtained estimates are 
also independent of e > 0. So, proceeding as in step 4, we can conclude the existence of 
{p,u, Pi, (32) and a subsequence of {p^^u^, PeUe, PeUeUe), still indexed by e, such that the 
following convergences hold true as e — )• 0: 

p, ^ p weakly - * in L°°(0, T; LP(0)), (2.35) 

p, ^ p strongly in C7([0, T]; iy-^'P(f7)), (2.36) 

u^-^u weakly in L^{<d, T; H^{n)), (2.37) 

{weakly in L'^{0,T; (Q)), 
weakly - * in L°°(0, T; L^hT (0)), (2.38) 
Strongly in L2(o,r;I^ ^'i^(J7)), 
PeUeUe (32 Weakly in ^^/^(O, T; l2(0)). (2.39) 

Moreover, using (j2.17p and the continuity of the product mapping from ff^(il) x ^¥"^'"^{0) 
to VF~^'2(f2), we obtain 

eVpe strongly in L^{Q,T]L'^{n)), 
eUeAp^ weakly in L^(0, T; l^-^'i (17)). 

Working as in the end of step 3, one can prove that /3i = pu and /32 = puu. Moreover, 
from above convergences, in analogous way as we did in step 5, we can pass to the limit 
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in the regularized problem (|2.3p . The only difference here is that the two terms involving 
the parameter e vanish. 

Thenceforth, we obtain that equation dtp + div {pu) = is verified in T>'{Q). From 
(|236]) we have that po^^ = Pe{0) p(0) in W-^'P{n). However, as po,e Po in LP{n), 
when 6 < p < oo, and />o,e ^ Po weakly-* in L°°{il), when p = oo, then p{0) = po at least 
in VF-^'P(fi). 

Finally, it remains to verify the energy inequality (12.2^ . As usual, integrating (I2.20p 
over (0,t), multiplying the result by <^ G 25(0, T), (p > 0, integrating over (0, T), and 
passing to the limit using convergences (I2.28p - (j2.33p . we arrive at the following inequality 



(hp,\u,Wt)\\i+ [ [2i^\\Du,is)f2 + 2ua [ \u,\^]ds] <P{t)dt 

Jo Jan J 

<\\ ^^dx / + [ ( [ I Pef-UA <P{t)dt. 

^ Jq Po,£ jo Jo \Jo Jn / 

I |2 I |2 

As e — )■ 0, we have ^^^^"^ dx — )• '^^dx < C. From convergences (|2.35p - (|2.39p . as e — ;> 
we find 



T /I 







2llPl^l'Wlli 



+ / [27^111)^(5)11^ + 21/0 / \u\^]ds) ct){t)dt 
Jo Jdn / 



|2 rT pT / pi 

2 Jo. Po 



< \ I ^^dx I 4>{t)dt + I [I I pf-u] 4>{t)dt, 
10 Jo \Jo Jn 



for any (/> G P(0,T), (p > 0. This yields the energy inequality (j2.2p which implies in 
particular that G (0, T; (J]) ) . □ 



Remark 2.2. Observe that, from the regularity of p and pu, equation dtp + div (pu) = 
holds in 

W-^'°°{0, T; LP{n)) n L°°(0, T; W~^'^{n)). 
Remark 2.3. Concerning the pressure, observe that 

pu E L'^{0,T;L^{n)) =^ dt{pu) G VF"^'~(0, T; (17)), 
puu G L^/2(0,r;L2(O)) ^ di?;(puu) G L'/^{0,T;W-^'^{n)), 
u G L2(0, T; F^(!^)) ^ An G L2(0, T; W'^'^in)), 
pf G Li(0,r;LFFT(f])) ^ p/ g L^O, T; W^-^'FFT(fi)). 

From i/ie De Rham theorem there exists a distribution vr G Vl^~^'°°(0, T; (J7)) suc/i 
that 

dt{pu) + div{puu) - uAu + Vvr = p/ in W-^'°°{0, T; W~^'^{n)). 
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3. Vanishing viscosity limit 

In this section we establish the convergence of a weak solution of the non-homogeneous 
Navier-Stokes equations with Navier boundary conditions to the strong solution of non- 
homogeneous Euler equations when the viscosity coefficient goes to zero. To this end, 
consider the following limiting problem 

dt{pu) + div {puu) + Vvr = pf in Q, 

div n = in Q, 

dtp + div (pu) =0 in Q, 

u ■ n = on S, 

/o(0) = po in n, 

n(0) = uq in il.. 



(3.1) 



We recall a result of existence of strong solutions to the nonhomogeneous Euler system 
(|3.ip . For another results on this subject see [TSj \2\ [2l]. 



Theorem 3.1. (see |34j ) Let Q with boundary dQ smooth enough, and p > 3. Assume 
that pq G VF^'^(il), < p^ < Po < p* < oo, for some constants p*,p*; moreover, assume 
that UQ G W^'P{n), uo-n = on dn, diviuQ) = 0, f € L^{0,T;W^'P{n)). Then there 
exists a time G (0, T] such that problem (j3.ip has a unique solution (p, u, vr) which 
satisfies 

a) p e C{[0,n];W^'P{n)) n C\{0,n];W^'P{n)), 0<p^<p<p* <oo, 

b) u e c{[o,n];W^'P{n)) nw^'^{o,n;w^^p{n)), 

c) ire L\0,n;W^^P{n)). 

Moreover, iff G C{[0,T]-W^^P{n)) thenue C^{[0,T,]-W^^P{n)) andn G C{[0,T,];W^^P{n)), 
hence {p, u, vr) is a classical solution. 

Remark 3.1. Observe that, since p > 3, by the Sobolev embedding we have that Vp G 
L~(Ox(0,T*)), Vu G L°°(Jlx(0,r*)), ut G L\0,T,; L°°{n)) andVn G L^{0,T,- L°°{n)). 

Let {u'^, p^) be a weak solution of 

' dtip^u") + div {p^u^u") - vl^u"" + V-n" = p" f in Q, 

div n'^ = in Q, 

dtp"" + div {p^'u'') = in Q, 

u'^ ■ n = on S, 

[Z)(n> + an1,,„ = on S, 

p-(o) = p'^ in n, 

{p''u''){0) =v^ inJl, 

with < Pq G LP(Q) given by Theorem 12.31 This weak solution satisfies 



(3.2) 



^-\\y^)u'^{t)\\l + 2:.aj^ 



an 



1 
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and 



< - 



u\\2 
2- 



(3.4) 



We now state the main result of this section. 



Theorem 3.2. Under the hypotheses of Theorems and \2.S[. Let {u,p) be the solution 
of ( fg.ip obtained in Theorem \3.1\ and let {u'^,p'^) be the one of \3. S^) given in Theorem 
2.3[ Assume further that < /3* < < p*, for all > 0, where and p* are the same 
constants of Theorem \3.1[ Then there exists C > independent of u > such that, for 
all t E (0,T*], the following inequality holds: 



+ Mt)-p' 



\\u{t) - U 

< C 

In particular, if 
L^{0,T^;L^{i})), then 



PqUo 



PqUq 



2 
2 

2 rTt 

+ WPO- Po\\2 + ^ 
2 Jo 



+ 



r 



I 2p 

p-1 



+ IIpo -Polli 



sup (\\u{s)-u'^{s)h + \\pis)-p'^{s)h) 

0<s<t 



as — 7- 0, and f^ converges to f in 
as i/ ^ 0. (3.5) 



Remark 3.2. We stress that (|3.5|) implies in particular that 

sup \\^/])u — \fp^u''\\2 — > as u ^ 0, 

0<s<t 

since we can write 

^pu - ^u'' = (Vp - \/p^)^i + ^[y- - u'"). 

Moreover, by interpolation, we conclude that, for p > 2, 



sup llp-p lip 

0<s<t 



as V ^ 0. 



Let us proceed with the proof of Theorem 13.21 
Proof. The differences u = u — u'^ , a = p — p^ , q = tt — tt^ satisfy 



p^'luji + u" ■ Vuj\ + Vq 
-p^co ■ Vu - (j{ut + u-Vu)- lyAu" + a/ + p^if - 

div oj = 
dtCF + u^ ■ Vfj = —uj ■ Vp 
OJ ■ n = 
a(0) =po-p'o 
w(0) =uo-u'^ 



in Q, 
in Q, 

in Q, (3.6) 
on S, 
in ri, 
in fi. 

Formahy, by multiplying third equation in (j3.6p by a and integrating in space and time 
we obtain 
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< 7Tlk(o)ll 



uj ■ V pa, 



Jn 



(3.7) 



where we have used that div{vy) = and u'^-n = on dO, to deduce that {u^ ■ Va) a = 
0. In fact, to justify the previous inequahty we proceed in this way: consider the energy 
inequahty (|3.4p and obtain other three inequahties: Ji by multiplying the continuity 
equation (jS.ip ci for p by p, J2 by multiplying the same equation (jS.ip ^ by p^ and J3 by 
multiplying the equation (|3.2p q for p'^ by p. One obtains (|3.7p by doing (j3.4p + Ji — J2 — ^3- 
From (|3.7p . by using the Holder and Young inequalities we arrive at 



kWlli<Ho)lli+ / llvpiioodlo^lli + lHli) 



(3. 



Formally, by multiplying the first equation in (j3.6p by w, and integrating in space and 
time, and finally by using integration by parts (Lemma 12. 2p in the Laplacian term, one 
gets 



1 

< - 
- 2 



i + 2ua I [ u" ■w + 2iy [ [ Diu") ■ D{w) 
Jdn Jo Jn 

t r ft 



V, 







p'^oj ■ Vu ■ UJ 



2 ./o Jn Jo Jn 

t r ft 



a{ut + u ■ Vu) ■ oj 



+ I / af-u:+ / / p-{f-n.u:. (3. 
10 Jn Jo Jn 



To justify the previous inequality, we first obtain three inequalities: Ii by multiplying the 
momentum equation ()3.ip i by u, I2 by multiplying the same equation (j3.ip i by u'^ and I3 
by multiplying the momentum equation (j3.2p i for u'^ by u. Next, we consider the energy 
inequality (j3.3p and do (j3.3p + Ii — I2 — I3 to arrive at (j3.9p . 
Proceeding as in [13], we rewrite 



2iya [ [ u" -w + lu [ [ Diu") ■ D{w) 
Jo Jdn Jo Jn 



2ua 



Jdn 



+ +2v 



Jn ^ 



\u\ 



Jan 



\D(u) 



Jn 
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Thus, by using Holder and Young inequalities, we have 



1 



Jo 



P^^Wl + / ll'Wt||oo||cr||2||w| 




JdQ. 



+ \\u- Vu||oo||cr||2||w||2 + 

Jo 

I ii/iiooikii2iio^ii2+ / w^hpWf - r\\2^\\^A2 

Jo Jo 



+ 



< 



PqUq 



vo 



+ / l|vu||oo||^/p^w||i 



. + ||n-Vn|U + ||/||oo)(lk||i + ||cT||2) 

+ Cv I + C||Vp^||ioo(o,r.;L2p(n)) / \\f - f\\^\\\^^h 

Jo Jo 

where C = C{p*,a, Q) is independent of u and we have used that / \u\'^ < C{Q)\\u\\ffi(^Qy 

Jan 



Notice also that 

1 1 / 771 1 II u w^/ 

WVP llL°°(0,T.;L°°(n)) — \\P \\l°°{0 

Adding (j3.8p and using ()3.10p . we obtain 



i/||l/2 ^ 

T.;L°°(Q)) — llf^Ulloo 



< C(p*)i/^ (3.10) 



< 



PqUo 



+ \Wm\l + Ci^ [ Mmin) + C [ ||/-ril^||a; 

Jo ^ ' Jo p-^ 

+ C [ (||Vp||oo + ||Vn||oo + ||ut||oo + ll-u • Vu||oo - 
Jo 

Now we apply a Gronwall type inequality [271 P-360] to obtain 



oo)(ll^lli + lkll2) 



Ht)-u^ml + Mt)-p''it)g 



< C 



P>o - 



+iiPo-poiii+^^ r 

Jo 



+ 



J 11^ 

p-i 



where C is a positive constant depending on 



p^,p*,a,n, ||Vp||ii(o,T,;L°°(n)), l|Vn||ii(o,T.;L°°(t7)), 
ll'"t||Li(0,T*;L°°(Q))i IKIlL°°(nx(0,T.)) and || /||Li(0,T,;L°°(n)) • 

Due to regularity of p, u and / in Theorem 13.11 (see Remark 13. ip . all the above norms are 
finite. The proof is then complete. □ 
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